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Abstract 



In the present paper, we propose a mean field approach for spin ladders 
based upon the Jordan- Wigner transformation along an elaborately ordered 
path. We show on the mean field level that ladders with even number legs 
open a energy gap in their low energy excitation with a magnitude close 
to the corresponding experimental values, whereas the low energy excitation 
of the odd-number-leg ladders are gapless. It supports the validity of our 
approach. We then calculate the gap size and the excitation spectra of 2-leg- 
ladder system. Our result is in good agreement with both the experimental 
data and the numerical results. 
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I. INTRODUCTION 



The study on low dimensional Heisenberg anti-ferromagnetic model is one of the most ac- 
tive research field in condensed matter physics .Haldane conjectured that for integer-spin 
one dimensional anti-Ferromagnetic chain a energy gap exists in the low energy excitation 
spectrum , but for half integer-spin case the excitation spectrum is gapless. The spin 1/2 
anti- Ferromagnetic chain can be solved exactly by Bathe ansatz The excitation spec- 
trum is found to be gapless. The measurement on the realistic ladder material such as 
S'rC'U203(two-leg), S'r2CM305 (three-leg), or(VO)2-P2C'7 (two-leg) |[TD| shows that the spin 
excitation gap is opened in the ladders with even numbers of legs, while for ladders with 
odd numbers of legs, no gap is found in spin excitation. This conclusion is predicted by 
the early numerical calculation [1^ and was explained qualitatively by Khveshchenko [§. 
In Khveshchenko's explanation a topological term appears in the effective Hamiltonian of 
the long wavelength dynamics in odd-leg ladder and is absent in even-leg ladder. Recently, 
G.Sierra |l| has mapped the ladder problem onto a effective one dimensional non-linear Sigma 
model. An extra topological term appears in odd-leg system and is absent in even-ladder 
system. This difference between the odd and even- leg ladders is essentially an extension for 
the difference between the half integer and integer Heisenberg chains. 

We also have the experience from the 1-D quantum Heisenberg spin chain, that it is 
not so trivial to incorporate the subtle physics of the topological term in a mean field 
approach. For two-leg ladder, the existence of energy gap in spin excitation for nonzero 
inter chain coupling J' is confirmed by various method such as Lanczos, quantum Monte 
Carlo |lT^ , p!5[| ,renormalization group 0, variational method 0, strong coupling expansion [0, 
spin liquid mean field approach , as well as the mean field approach based on the Jordan- 
Wigner transformation |T^. And for J=J' (J is the exchange coupling along the chains) 
which is the physical parameter of the real ladder compounds, the energy gap obtained by 
numerical calculation is 0.5 J which is very close to the experimental value. The spin wave 
excitation spectra has also been obtained by the numerical calculation which shows the 
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minimum of the spectra is located at the wave number vr [p4l|T5i. This was also confirmed 



by the recent neutron scattering experiment 11 



A mean field treatment of the two-leg Heisenberg ladder with the application of Jordan- 
Wigner(J-W) transformation has been proposed by M.Azzouz and et al [^]. In their ap- 
proach the Jordan- Wigner transformation is introduced to map the spin 1/2 system to the 
spinless Fermion system. The gap is obtained in their mean field approach which is about 
0.7J' under the case of J=J' which does not fit well with the experimental value. The exci- 
tation spectra is also calculated in their approach, which contents a minimum at the wave 
number tt. But the shape of the spectra is not consistent with the numerical result which 
predicts the maximum of the spectra locating between the and 7r/2 [p!5[] . 

In the present paper, we propose a mean field approach also based on the J-W trans- 
formation ||l^. Our mean field approach is quite different from the approach used in paper 
||T6| . When performing J-W transformation in the ladder system, one must put the sites in a 
queue. Then the spin operator can be expressed as = c^e*'^**-*\ where is nothing but 
the summation of number of the sites from — oo to the (i-l)th site in the particularly ordered 
queue which are occupied by spinless Fermions. The difference of particularly elaborated 
queues used in our approach and that in paper are shown in Fig.l.(a) and (b). In our 



approach the sites in the odd number rungs are labeled from upward to downward but from 



downward to upward in the even number rungs(shown as Fig. 1(a)), whereas in paper |T6 
all rungs are labeled from upward to downward, (shown as Fig. 2(b)). 

We perform the J-W transformation with such a specially ordered queue. The advantages 
of our approach are as the following. First of all if we replace the phase of the hopping term 
by its average value, we can easily obtain the mean field Hamiltonian for n-leg ladders in 
which the effective hopping terms of the spinless Fermions have alternative signs along both 
the chain and rung directions. This mean field Hamiltonian is similar with the one used in 
paper |T^ for 2-leg case. But in their approach a further assumption, i.e. adding a vr-fiux 
to each plaquette is needed to obtain such kind of mean field Hamiltonian, whereas in our 
approach the mean field Hamiltonian is obtained directly by replacing the phase factor in 



the hopping terms with their average values. We can then easily show that the spin gap is 
only opened in even-leg ladders. In certain degree, it supports the validity of our way of 
Jordan- Wigner transformation construction. 

Secondly we can treat the phase factor in the hopping terms more carefully by introduce 
a self consistent procedure for 2-leg ladder. We find that both the gap magnitude and the 
spin excitation spectra are in good agreement with the numerical and experimental results in 
the case of J'=J which is of the physical parameter for the real ladder compounds. Actually 
he energy gap is found to be 0.46J in the case J'=J which is very close to the experimental 
result 0.47 ±0.2 J. The spin wave excitation spectra obtained by our mean field calculation 
is also consistent with the numerical result with the minimum at the wave number vr and 



the maximum at the wave number 0.3567r. Compared with the approach used in [|T^], our 
mean field approach is much better in the case of J' > 0.5 J. In the weak coupling regime 
our approach works not so good, the energy gap is even persisted in the case of J' = 0. In 
spite of the unsatisfactory aspect in the weak coupling regime, our approach is still valuable 
because it works very well in the intermediate and strong coupling regime, which has the 
experimental correspondence. 

In section II, we calculate the spin excitation gap of the spin ladders with various numbers 
of legs and show that in the mean field level the energy gap only exists in ladder system 
with even number legs. In section III, we propose a more carefully treatment of the 2-legs 
spin ladder. The spin excitation gap, as well as the excitation spectra is calculated. Finally, 
we make the conclusion remarks in Sec. IV. 
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II. Mean field treatment for the spin gap in n-leg ladders 

We begin with the 2M-leg anti-ferromagnetic Heisenberg ladder Hamiltonian: 



2M-1 ^ ^ 2M ^ ^ 

H = J' ^ ^ Si^p ■ Si^p^i + J ^ ^ Si^p ■ Si^i^2M+i-p 

i p=l i p=l 

2M~l 2M 

= ^ X/ X/ ^i,p ' ^i,p+l + "^X/ X] ^i,p ' '^i+l,2Af+l-p+ 
i p=l i p=l 

Jl 2M-1 J 2M 

^ X! "^tp ' "^i^P+l + "2 51 X! "^tp ' ^i+l,2M+l-p + H.C. (1) 

i p=l i p=l 

In the above Hamiltonian, 'i' represents the site position along the chains and 'p' represents 
the 2M sites of different chains coupled by the inter-chain coupling constant J'. As shown in 
Fig.l, p is labeled from upward to downward at the even sites while downward to upward 
at the odd sites. This is different from what used in the paper of Azzouz et.al. |16|. Then 
we introduce the generalized J-W transformation: 

^p^-^p > ^'=1 ^1=1 (2) 

in which c is the spinless Fermion operator. The summation in the phase factor is nothing 
but the number of the occupied sites before the ith site along the particular queue shown in 
Fig. 1(a). Then the quantum spin 1/2 Hamiltonian can be mapped onto a spinless Fermion 
Hamiltonian as: 



2AI-1 2 I 2M ^ 2 

H = J' ^ {1^- ctp(^i,p) ■ {-j^- ctp+l(^i,P+l) + J ^{2~(4,p(^i,p)'{2~(^t+l,2M+l-p(^i+l,'2M+'i--p) + 
i,p=l *)P=1 

T/ 2M+1 T 2M 

y E «pC.P+i + H.C.) + E«pC,+,,2M+i-pe-^*(^) + H.C) (3) 

i,p=l i,p 

where 

2m 

i=p+l 

In our mean field approach, we replace rn^i by < ni^i >. For the present study, we further 
assume that the finite magnetization in this system is not possible because of the strong 



quantum fluctuation. This is reasonable for systems with its leg number much less then the 
site number of each chain. Then we have < S^i >=< 1/2 — cfiCi^i >— which implies < 
n^i >=0.5. Consequently, the phase factor in eq.(3) can be approximated by:$(p) = 7r(2M — 
p) . Moreover we decouple the four fermion interaction term in the above Hamiltonian by the 
Hatree-Fock approximation. Finally the mean field Hamiltonian of the spinless Fermions 
has the expression: 



J' 



-I 2M+1 



2M 



^ = 2" ^ ^<v''^'^+' + + E«pC.+i,2M+i-p(-l)^+^ + H.C.) 

i,p=l i,p 

If we introduce a Fourier transformation for the site indices, we have: 

2M+1 J 2M 

H = a J2 i4,pCk,p+i + H.C.) + -^(i7fc(-l)P4pCfc+i,2M+i-p + H.C.) 



(4) 



(5) 



k,p=l 



k,p 



where a — ^ Ik — —Jsm(k). The Hamiltonian then can be written in a form H 
Efe C+M^)Cfe with h(k): 



h{k) 



a 

a a .. 

a a .. 

.. .. a a 



?7 



— Z7 



Z7 



^7 a a 

.. —i'j a a 

^7 a a 

y —ij a 

The above matrix contents 2M eigenvalues for a given wave number k corresponding to 
the 2M individual bands separated by gaps. It can be proved straightforwardly that for a 



/ 
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given wave number k , half of the eigenvalues are less than zero and other half of them are 
great than zero. Furthermore we can also prove that zero is not an eigenvalue of the above 
matrix for any non zero J'. We will prove the two statements in the Appendix. This result 
shows clearly that half of the energy bands of the spinless Fermions is below zero energy 
and another half is above zero. The energy gap between them is nonzero because zero is not 
the eigenvalue of the above matrix for nonzero J' and arbitrary k. Assuming there are no 
self-magnetization in one dimensional system, only the lower half of the states is occupied 
by spinless Fermions in ground state. Therefor such spinless Fermion system is very similar 
with the traditional insulator in which the valence band is fully occupied and the conductive 
band is fully empty in ground state. So a nonzero minimum energy is needed to excite 
the system from the ground state which indicates a spin excitation energy gap. For spin 
ladders with odd number legs there exist odd numbers of energy band. Since only half of the 
states is occupied in the ground state there must exist at least one band which is partially 
occupied. This picture is very similar with the traditional conductor in which exists at least 
one partial occupied band (conduction band). Then the low energy excitations are gapless 
for odd-number-leg ladders. 

Moreover, we calculate the gap size of the 2,4, 6, 8, and 10-leg spin ladder, the results are 
shown in Fig. 2 together with the experimental value. Although our approach is quite rough. 



the result is in good agreement with the experimental value ||T0|,|Tl],|T2|,|T3[ . 

The spin excitation spectra can also be obtained from the above mean field approach, the 



spin wave dispersion is y ( y)^ + sin A;^ for the two-leg case. It has two energy minimum, 
one is at and another at vr. The shape of the spectra is not consistent with the numerical 
result which has the spectra minimum located only at wave number vr and the maximum is 
near O.SSGvr. This is due to that the treatment for the phase term of eq.(3 ) is too rough. 
In the next section we will propose a more careful treatment of the phase factor in eq.(3). 
We can then obtain a much more improved spin excitation spectra, which is very close to 
the numerical results. 
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II. The mean field theory of 2-leg ladder 



For the two-leg case, we can introduce two bipartite lattice labeled as a and (3. Following 
the same procedure shown in the above section, the Hamiltonian for spinless Fermions 
becomes: 



i i i 

^ E(«^^A + h.c.) + ^ J2{at/3,^,e^-^^tPi+<.'^^^^^ + h.c.) + ^ + h-c) (6) 

it i 

In our mean field approach, different with the simple treatment used in the previous section, 
we first replace the phase factor in (6) by it's average value: 



= -4xiX2 

where we define: 

Xi =< Pt(^i+i > X2 =< a+^iPi > Xo =< A^ttj > 

Then the Fermion-Fermion interacting term (1/2 — Q;^Q;j)(l/2 — jS^iPi+i) can be factorized 
as: 

(1/2 - a+a,)(l/2 - /^^lA+i) = \- XoP^cx, - x^a+A + Xo Xo 

We decouple the other two interacting term in the same manner and obtain the following 
mean field Hamiltonian of spinless Fermions. 

Hmf = Y^7kOit(3k + h.c. (7) 

k 



where 

2 ^2 



7fe = ( - J'Xo) + d- 2J|x2r - Jxi - JX2) cos(A;) ) + 
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ism{k) (^Jx2 - Jxi - 2^|X2p - ^ 
Then the above Hamiltonian can be diagonahzyed as: 

Hmf = E Ek {at^k - ~Pth) (8) 

k 

in which : 

Ek = 

i ((y - ^'^o) + - 2J\X2? - Jxi - JX2) cos(A;)) + sin2(A;) i^Jx2 - Jxi ~ 2J\x2? - ^ 

The three parameters X15X2 and xo is determined selfconsistantly. The gap size obtained 
within the present approach is shown in Fig. 3 compared with the numerical results. Our 
results fits quite well to the numerical results in the parameter regime J' /J > 0.5. For 
the case of J' = J, the three parameters are found to be Xi = — 0.188J X2 = 0.237J and 
Xo = 0.3867 J, and the gap is found to be 0.46J which is very close to the experimental value 
(0.47 ± 0.2) J [|l|]. In strong coupling limit (J' >> J) our result fit well with the result 
obtained by strong coupling expansion which shows A/J'^1 when J'/ J —>■ 00. But in 
the regime J'/ J < 0.5 our results deviated from the numerical results, and a nonzero gap 
persists even at the case J'=0. So our mean field approach is valid only in the intermidate 
and strong coupling regime. In the weak coupling regime our mean field picture breaks due 
to the over estimation of the inter-chain interaction. Since the phase factor in the hopping 
term in equation(6) is replaced by its average value, it makes the hopping term within 
one chain being strongly modified by the motion of spinless Fermions in the other chain, 
which is not valid for the weak coupling regime. We believe this approach is valid when the 
inter-chain coupling is in the order of unit, but is not valid for the weak coupling case. 

Another advantage of the present mean field approach is that in the case of J' = J 
it gives the same spin wave dispersion predicted by the numerical calculation as shown in 
Fig. 4. The minimum of the spectra is at the wave number vr, and the maximum is at the 
wave number 0.3567r. This result is in good agreement with the numerical result which has 
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a minimum at tt and maximum at O.Svr. We can also calculate the two-magnon continuum 
from our mean field theory. The two magnon continuum is proportional to 



Jdt< S'{q,t)S'{~q,0) > e-'^' 



which can be transformed into the density- density correlation of the spinless Fermions by a 
Jordan- Wigner transformation. Then the two magnon excitation can be viewed as particle- 
hole excitation of the spinless Fermions. The spectra of the two magnon excitation with 
several specific q number is shown in Fig.5. And the bottom(top) of the two magnon con- 
tinuum is just the minimum(maximum) energy of hole- particle excitation of the spinless 
Fermions for a given wave number. The result is shown by dashed line(bottom) and dot- 
ted line(top) in Fig. 4 which fits the numerical result quite well. 0,0. Compared to the 



numerical method such as DMRG [0, quantum Monte Carlo and Lanczos method |TJ,|15 
our mean field theory based on the Jordan- Wigner Transformation gives a more transparent 
understanding of the gap formation in even-number-leg spin ladders and the low energy spin 
excitation. 

The spin susceptibility is also obtained by introduce a magnetic field in the original 
spinless Fermion Hamiltonian, this term acts like the chemical potential: 

k k 

The magnetization m then has the expression as: 

^ = j^ < ^t<^i > - < A^A >) 

i 

And the spin susceptibility could be derived as: 

dm 

We calculate the spin susceptibihty in the case of J' = J in a wide range of temperature, 
the result is shown in Fig.6. Our result explains the temperature behavior of the spin 
susceptibility quite well, and it is again in good agreement with the numerical results which 
gives a maximum at T = 0.8 J. 

10 



IV. The Concluding Remcirks 



In this paper we propose a mean field approach for spin ladders based on the Jordan- 
Wigner transformation along an elaborately chosen path defined in this paper. We show 
that in the mean field level that spin gap is opened only in the even-number-leg-ladders 
and vanishes in the odd- number-leg-ladders. It gives a very simple picture of the formation 
and vanishing of the spin gap in the above mentioned two type of spin ladders. The spin 
ladders with even number legs formed a 'insulator' like band for spinless Fermions, whereas 
in odd-number-leg ladders the band structure of the spinless Fermions is 'metal' like. 

Then we take a more careful study of the 2-leg ladder. Particular for the J = J' 
case the magnitude of the gap found in our approach is in good agreement with both the 
numerical result and the experimental result. Further the spin excitation spectra and the 
uniform susceptibihty are also calculated based on our mean field treatment. The dispersion 
relation of the spin excitation spectra obtained by our mean field theory is very similar 
with the numerical result, which has its maximum locating between and tt. The uniform 
susceptibility is also consistent with the numerical results, which predicts a maximum at 
T = 0.8 J. 
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Appendix 

In the Appendix we prove that matrix h(k) in section II has the two following proper- 
ties: (i) If A is a eigenvalue of the matrix, —A is also a eigenvalue of it. (ii) Zero can not be 
a eigenvalue of the matrix with any nonzero J'. 

First we divide the Hermite matrix h(k) into its real and imaginary part h{k) = A + iB, 
in which the matrix A,B read as: 

Aij = aSij+i + aSij-i Bij = (-l)*+^(5j,2M+i-j7 

One can easily find that the matrix A and B satisfies some relations: 
(7) K ■ h{k) ■ K = -h{k) 

{II) = A B^' = -B B-^ = 7-25 B-^AB = -A 

in which Kij — Sij{—iy and KuKij — 5ij 

Based on the above equations we can prove the two properties straightforwardly. For (i)if 
we have h{k)x — Xx, we can multiply the matrix K to the both side of the above equation: 
K ■ h(k) -K-Kx^ XKx, then we have h{k){Kx) = -X{Kx). 

For property (ii), we can prove it as follows. First if B=0 the conclusion is obviously 
true because the determinant of matrix A is nonzero for J' ^ and this makes the equation 
Ax — can not be satisfied unless x=0. Next when B ^ we have : 

{A + iB)x — Q or equivalently B~^Ax — —ix 

For matrix C = B~^A, we have 

C+ = [B-^A]+ = ^-''A+B+ = -j-^AB = -{B-^){^-^B)AB = iB'^)A = C 

Then the matrix C is Hermite, and it can not has a imaginary eigenvalue, so the equation 
{A + iB) ■ X — Q can not be satisfied for nonzero x. 

In the above proof we used the relation (11). Then from the above paragraph we prove 
the two properties used in section II. 
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FIGURE CAPTIONS 



Fig. 1: To perform the Jordan- Wigner transformation the sites are put in a particular 
queue in our present study(a) and in paper |[T6[|(b). 

Fig. 2: The spin gap for n-leg spin ladders calculated in our mean field approach which 



is compared by the experimental value. |Tl],IT3 



Fig. 3: The solid line shows the spin gap obtained by the self consistent procedure for 
2-leg ladder as a function of inter chain coupling J'. The dashed line is the result of strong 



coupling expantion ||^. The squares shows the numerical result of paper |T^. The inset 
shows the results with J'/J less than 0.8. 

Fig. 4: The solid line is the dispersion of the spinless Fermions for 2-leg ladder calculated 
in our mean field approach. The dashed(dotted) line is the bottom(top) of the two magnon 
continuum. 

Fig. 5: The spectra of two-magnon excitation with q = O.l7r(dotted line), 0.57r(solid 
line), TT (dashed line). 

Fig. 6: The spin susceptibility of 2-leg ladder in the case of J = J'. 
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Fig.3 
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